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Abstract. We give an overview of circumstellar interaction in young Type II supernovae, as seen through the eyes of very- 
long-baseline interferometry (VLBI) observations. The resolution attained by such observations ( < 1 mas) is a powerful tool 
to probe the interaction that takes place after a supernova goes off. The direct imaging of a supernova permits, in principle, to 
estimate the deceleration of its expansion, and to obtain information on the ejecta and circumstellar density profiles, as well as 
estimates of the magnetic field intensity and relativistic particle energy density in the supernova. Unfortunately, only a handful 
of radio supernovae are close and bright enough as to permit their study with VLBI. We present results from our high-resolution 
observations of the nearby Type II radio supernovae SN 1986J and SN 2001gd. 



1. Introduction 

Radio supernovae are mostly associated with core-collapse su- 
pernovae (Type II and Ib/c). Here, we will discuss only Type 
II radio supernovae, although most results apply as well to 
Type Ib/c radio supernovae. The progenitors of Type II su- 
pernovae are surrounded by a high-density wind, n ~ 3 x 
10 7 M_5 ris v^ 1 cm -3 , where Ms is the mass loss rate of the 
presupernova wind in units of 10~ 5 M yr~', vi is the pre- 
supernova wind velocity in units of lOkms , and r^ is 
the supernova shock radius in units of 10 15 cm. In the stan- 
dard interaction model for radio supernovae (Chevalier 1982 
Nadyozhin 1985 1, a forward shock at velocities thar largely 
exceed ~ 10000 km s _1 propagates into the circumstellar 
medium (CSM) when the blastwave reaches the dense, ionized, 
slowly expanding wind. In addition, a reverse shock moves 
back into the stellar envelope at speeds of ~ 1000 km s _1 rela- 
tive to the expanding ejecta. As a result, a high-energy-density 
shell is formed. 

The brightness temperature inferred from VLBI images of 
SNe is very large, indicating that the emission is of non-thermal 
origin, namely synchrotron radiation from relativistic electrons 
(Ne = No E~ y ) in the high-energy-density shell. As it expands, 
a radio supernova quickly increases its radio brightness with 
time, due to the increasingly smaller electron column density 
in the line of sight. When the optical depth at cm-wavelengths 
has dropped to about unity, the supernova reaches its maximum 
of emission, after which the emission monotonically decreases. 
If synchrotron self-absorption effects are negligible, then the 
only relevant source of absorption is free-free absorption by 
thermal electrons in the CSM, and the radio luminosity evolves 
as L Y oc R 2 AR No B 1 ^ v _z r e _T ^ (r ffjV oc v" 21 ). Thus, a radio 
supernova is first seen at high frequencies and, as times goes, 
it becomes visible at increasingly lower frequencies. We see 
that the flux density monitoring of radio supernovae may yield 
important information on such fundamental physical parame- 
ters as magnetic field intensity and relativistic particle energy 



density in the supernova shell, as well as information on the 
relevant absorption and cooling mechanisms. Since the radio 
peaks of SNe are usually of the order of a few milliJansky, only 
bright, nearby radio supernovae are good VLBI targets. 

The density of the supernova ejecta can be approximated by 
a steep power-law profile (p e j oc r~"). If the mass loss param- 
eters (M, v w ) stay approximately constant up to the explosion, 
the CSM density profile is then p csm oc r~ 2 . The dense shell 
formed at shock breakout decelerates as it interacts with the 
CSM, and it can be shown that the expansion of the shell radius 
follows a power-law with time: R oc f"; m — (n - 3)/(n - 2). 
More generally, p csm oc r~ s , leading to m = (n - 3)/(n - s). 
Therefore, by carrying several VLBI observations of the same 
supernova, we may be able to infer the value of its deceleration 
index, m, and indirectly the values of n and s that characterize 
the ejecta and CSM density profiles. 

2. SN 1986J in NGC 891 

SN 1986J in NGC 891 is one of the most radio luminous SNe 
ever discovered. The precise date of its explosion is not known, 
but on the basis of the available radio and optical data SN 1986 J 
was estimated to have exploded around the end of 1982, or 
the beginning of 1983 (e.g., Weiler, Panagia & Sramek 1990). 
Based upon its large radio luminosity, Weiler et al. 1990 sug- 
gested that the progenitor star was probably a red giant with 
a main-sequence mass of (20 - 30) M that had lost material 
rapidly (M ^ 2 x 10~ 4 M yr -1 ) in a dense stellar wind. VLBI 
observations made at 8.4 GHz at the end of 1988 by Battel 
et al. ( 119911 hereafter B91) showed that the radio structure of 
SN1986J had the form of a shell, in agreement with expecta- 
tions from the standard interaction model, with the minimum 
of emission located approximately at its center. Those authors 
claimed the existence of several protrusions at distances of 
twice the shell radius, and with apparent expansion velocities as 
high as ~ 15000 km s _1 . Since these protrusions were twice as 
far as the mean radius of the shell, it then follows that the main 
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Fig. 1. Left: 5 GHz VLA map of the galaxy NGC 891 and its supernova SN 1986J made on February 21, 1999. The contours are 
(1, 2,4,8, ...)x 106 /iJy beam -1 . The peak of brightness of the map corresponds to the supernova and is ~ 8.2 mJy beam -1 . Right: 
Global VLBI map of SN 1986J on February 21, 1999. The contours are (3,5,7,9,1 1,13,15,17,19) x 56 fiJy beam 1 , the off-source 
rms. The peak of brightness of the map is 1.13 mJy beam -1 and the restoring beam (bottom left in the image) is 1.3 x 0.9 mas 2 
at a position angle of -13°4. In both panels, north is up and east is left. 



bulk of the shell expanded at roughly 7500 km s . Such pro- 
trusions have been successfully invoked by Chugai ( 1993) to 
explain the coexistence of velocities smaller than 1000 km s -1 
implied from the observed narrow optical lines (Leibundgut 
et al. 11991b . and the large velocities indicated from the VLBI 
measurements. 

We used archival VLA and global VLBI observations of su- 
pernova 1986J at 5 GHz, taken about 16 yr after the explosion, 
to obtain the images shown in Fig. [0 The right panel corre- 
sponds to the 5 GHz VLBI image of SN1986J. It shows a highly 
distorted shell of radio emission, indicative of a strong defor- 
mation of the shock front. The apparent anisotropic brightness 
distribution is very suggestive of the forward shock colliding 
with a clumpy, or filamentary wind. Note that there are several 
"protrusions" outside the shell, though just above three times 
the noise level and at different position angles from those previ- 
ously reported by B91. Therefore, these protrusions could not 
be real, but must be just artifacts of the image reconstruction 
procedure. If this is the case, the disappearance of the protru- 



sions seen in the previous VLBI observations (B91) would im- 
ply a change in the density profile of the circumstellar wind. 

At a distance of 9.6 Mpc (Tully 1998), 1 mas corresponds 
to a linear size of ~ 1.4xl0 17 cm « 0.05 pc. Based on 8.4 GHz 
VLBI observations, B91 found an angular size of ~3.7 mas for 
the shell of SN 1986J, 5.74 yr after its explosion (assuming it 
took place on 1983.0). The corresponding mean linear velocity 
would then be ~ 14700 km s -1 for the first 5.74 yr. However, this 
velocity applies only to the protrusions found by B9 1 , not to the 
shell. Indeed, the velocities reported in B91 were calculated 
for the protrusions, and assuming that these originated in the 
centre at the time of the explosion. These authors also pointed 
out that the protrusions extended from the centre to twice the 
radius of the shell, i.e., the protrusions were outside the shell. 
Therefore, a value of ~ 1 .85 mas for the angular size of the shell 
of SN 1986J at epoch 1988.74 is indicated, and a mean linear 
velocity of the radio shell of ~7400kms -1 is more appropriate 
to characterize the first 5.74 yr of the expansion of SN 1986J, 
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Fig. 2. Left: 8.4 GHz VLA map of the galaxy NGC 5033 and its supernova SN2001gd, from observations on 8 April 2003. 
The lowest contour is at 48//Jy beam 1 , corresponding to three times the off-source rms. The peak of brightness of the map 
corresponds to the nucleus of NGC 5033 and is 2.44 mJy beam 1 . The supernova is the bright point-like source northwards of 
the nucleus of NGC 5033. Right: (Top) VLBI map of SN 2001gd on 26 June 2002. Contours are drawn at 10, 20, 30, 40, 50, 60, 
70, 80, and 90% of the peak of brightness of 3.64 mJy beam 1 . (Bottom) VLBI map of SN 200 lgd on 8 April 2003. Contours 
are drawn at 10, 20, 30, 40, 50, 60, 70, 80, and 90% of the peak of brightness of 947/vJy beam -1 . In all panels, north is up and 
east is left. 



as has been previously noticed by Chevalier ( 1998 1 and Houck 
et al. dl998l . 

The angular size of the shell of SN 1986J on 21 February 
1999 is ~4.7 mas, equivalent to « 6.8 x 10 17 cm. Combining 
this angular size measurement with that obtained by B91 for 
epoch 1988.74 (9 « 1.85 mas), we obtain a mean angular 
expansion velocity of the shell between 29 September 1988 
(1988.74) and 21 February 1999 (1999.14) of « 0.14 mas yr 1 , 
which corresponds to a linear velocity of ~ 6300 km s _1 . If we 
assume that SN 1986J freely expanded for the first 5.74 yr of its 
life (r oc f,m = 1) and then started to decelerate, the expan- 
sion between the two epochs of VLBI observations is character- 
ized by m = 0.90 + 0.06. This mild deceleration contrasts with 
the case of several other supernovae, for which a strong decel- 
eration has been measured (SN1979C: Marcaide et al. 2002 
SN1987A: Staveley-Smifh et al. [T9951 Gaensler et al. 119971 
SN1993J: Marcaide et al. [T997l Barrel et al. 120001 . 

For a standard presupernova wind velocity, v H ,=10 km s _1 , 
the linear size of SN 1986J at epoch f=16.14 yr implies that 



we are sampling the progenitor wind about 11000 yr prior to 
its explosion. Since the mass loss rate of SN 1986 J seems to 
have been t 2 x 10~ 4 M Q yr _1 (Weiler, Panagia & Sramek 
1990), the swept-up mass must have been M sw » 2.2 M for 
a standard density profile of the progenitor wind (p cs oc r~ 2 ). 
Since the expansion of the supernova has not decelerated sig- 
nificantly between the two VLBI observations, it follows that 
the swept-up mass by the shock front must be much less than 
the mass of the ejected hydrogen-rich envelope, M env , as oth- 
erwise we should have observed a much stronger deceleration. 
In fact, momentum conservation implies that M en v ^ 12 M , 
significantly larger than M sw . If the hydrogen-rich mass enve- 
lope was as high as 12 M , this is a hint that the progenitor of 
SN 1986 J was probably a single, massive Red Super Giant (as 
previously suggested by Weiler et al. 1990), which lost mass 
rapidly, but managed to keep intact most of its hydrogen-rich 
envelope by the time of explosion. This result contrasts with 
the cases of SN 1993 J and SN 1979C, whose hydrogen-rich en- 
velopes had masses of 0.2-0.4 M (Woosley et al. 1994, Houck 
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& Fransson 1996) and ^ 0.9 M (Marcaide et al. 2002), and 
whose progenitor stars were part of binary systems. 

Since the radio emission is of synchrotron origin, we can 
estimate a minimum total energy (in magnetic fields, elec- 
trons, and heavy particles) and a minimum magnetic field for 
SN 1986J. If we assume equipartition, then the minimum total 
energy is (Pacholczyk 1970) = c 13 (l+k) 4/1 3/7 R 9/1 L 4/1 , 
where L is the radio luminosity of the source, R is a charac- 
teristic size, cn is a slowly-dependent function of the spectral 
index, a, <p is the fraction of the supernova's volume occupied 
by the magnetic field and by the relativistic particles (filling 
factor), and Ic is the ratio of the (total) heavy particle energy 
to the electron energy. This ratio depends on the mechanism 
that generates the relativistic electrons, ranging from A: a: 1 up 
to k — m p /m e m 2000, where m p and m e are the proton and 
electron mass, respectively. 

Based on snapshot VLA observations of SN 1986J carried 
out on 13 June 1999, we determined an spectral index of a = 
-0.69 + 0.06 (S v oc V) between 1.6 and 8.5 GHz. This value 
is very close to a — -0.67 obtained by B91, and is a typical 
value for supernovae that are in the optically thin radio regime. 
With a = -0.69 and S4.9GHZ = 7.2 mJy from our observations, 
we obtain a radio luminosity L v » 1.38 X 10 37 (D/9.6Mpc) 2 
ergs 1 for the frequency range 10 7 — 10 10 Hz. The value of the 
function co is approximately 2.39 x 10 4 . As the characteristic 
size for SN 1986 J, we take half the largest diameter of the shell, 
9 = 2.35 mas, which corresponds to a linear size ofR-D-9 a 
3.4 x 10 17 cm. With these values, the minimum total energy is 
then 

1.13 xlO-d+^^^D^ erg 

where 0=0.6600.66, 9-2. 3502.35 m as , and D=9.6Dg ^ Mpc. The 
value of the magnetic field that yields E^ m) is then equal to 

B min * 10.6 (1 + k) 211 <t>-% 9 2 f 5 D 9 fmG 

Since 1 < k < 2000, £ (tot) can have values between ~ 2 x 10 48 

mm 

and ~ 9 x 10 erg, while the corresponding values of the mag- 
netic field can lie between ~ 13 and ~ 93 mG (for (f> = 0.66 and 
D = 9.6Mpc). These values of the magnetic field for SN 1986J 
are in agreement with, e.g., those obtained for SN 1993J at sim- 
ilar radii. For example, Perez-Torres, Alberdi, and Marcaide 
l200T1 showed that B * 0.58 (r/3 x 10 16 cm) 1 G, which results 
in a value of ~ 50 mG for a radius of 3.4 x 10 17 cm. Since it is 
very unlikely that the magnetic field energy density in the wind 
is larger than its kinetic energy density, i.e., B 2 /8n < pv^/2, 
we can then obtain an upper limit for the magnetic field in the 
circumstellar wind of SN 1986J: 

B cs < (M v H ,) 1/2 r 1 ^ 0.78 (M-4 vio) 1/2 r\) mG 

where M=10" 4 M_4 M yr -1 , v w =10vio kms -1 , r=10 17 ri 7 cm, 
and we have assumed a standard wind density profile. For 
SN1986J, M = 2 x 10~ 4 N^yr 1 and r = 3.4 x 10 17 cm, 
and we obtain B csm 5 0.32 mG. These arguments suggest that 
the magnetic field in the shell of SN 1986 J is in the range 13-93 
mG, or about 40 to 300 times the magnetic field in the circum- 
stellar wind. Since a strong shock yields a fourfold increase 



in the particle density, the post-shock magnetic field is also 
fourfold increased. Hence, compression alone of the circum- 
stellar wind magnetic field cannot account for the large mag- 
netic fields existing in SN 1986 J, and other field amplification 
mechanisms are needed to be invoked, e.g., turbulent amplifi- 
cation (Chevalier 1982; Chevalier & Blondin 1995). The same 
conclusion was reached for the case of SN 1993 J (Fransson 
& Bjornsson 1998, Perez-Torres et al. 2001), where magnetic 
field amplification factors ~ 100 were found to be necessary. 

3. SN 2001 gd in NGC 5033 

NGC 2001gd in NGC 5033 was discovered by Nakano et 
al. J200U on 24.820 November 2001, although its explosion 
date is uncertain. The supernova had a visual magnitude of 
14.5, and is located ~3" north-northwest of the nucleus of 
NGC 5033. Nakano et al. (2001 1 reported the following posi- 
tion for SN 2001gd: a=13hl3m23s.89, 5=+36°38' 17".7. They 
also reported that there was no star visible at the above position 
on earlier frames taken from 1996 to April 2001. A spectrum 
obtained by P. Berlind ten days later (Matheson et al. 2001 1, 
showed SN 200 lgd to be a Type lib supernova well past max- 
imum light. Matheson et al. (2001 ) pointed out that the spec- 
trum was almost identical to one of SN 1993 J obtained on day 
93 after explosion Matheson et al. ( I2000> . Since SN 1993J was 
a strong radio emitter, it was natural to expect SN 200 lgd to be 
also strong in radio. Stockdale et al. l|2D02 1 detected SN 2001gd 
on 8 February 2002 at cm-wavelengths with the Very Large 
Array (VLA), thus confirming the suggestion that SN 200 lgd 
should be a strong radio emitter. Their continuous monitoring 
of SN 200 lgd since its first radio detection has confirmed that 
SN 200 lgd is very similar in its radio properties to SN 1993 J 
(Stockdale et al. 2003 1, as suggested by their optical similarity. 

We observed SN 200 lgd on 26 June 2002 and 8 April 2003 
at a frequency of 8.4 GHz, using an earth-wide global VLBI ar- 
ray (Perez-Torres et al. 2004). Unfortunately, we were unable 
to resolve the fine radio structure of SN 2001gd (Fig.|2j, despite 
the milliarcsecond resolution of our VLBI array. Therefore, 
we had to resort to non-imaging data analysis to make fits to 
the visibility data and estimate the true supernova size at each 
epoch. We assumed spherical symmetry for our model fits, and 
used the following models: (i) an optically thin sphere; (ii) an 
uniformly bright, circular disk; and (iii) an optically thin shell 
of width 25% the inner radius. Formally, all three models gave 
equally good fits to the (u, v) data, so in principle we could 
not rule out any of them. Although we suggest the uniformly 
bright, optically thick disk might be ruled out as a physical 
model in 2003.27, since the supernova was optically thin at 
8.4 GHz, we have left it for comparison purposes. The system- 
atic increase in angular diameter in all three models indicates 
that the supernova is expanding. 

At a distance of 13 Mpc, lOOjuas correspond to a linear 
distance of 2xl0 16 cm. The average expansion velocity of the 
supernova can then be written as 22600Z)i3 9r l km s _1 , where 
D=13 Do Mpc is the apparent distance to the supernova, 9 the 
angular radius in //as, and t the time since explosion, in days. 
Hence, the average expansion speed of SN 200 lgd on 26 June 
2002 (t = 296 d ), ranges from v exp = 12700+2600km s" 1 
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SN 1979C 1 SN 1986J 2 SN2001gd 3 SN 1993J 4 



Distance (Mpc) 


16.1 


9.6 


21.6 


3.63 


Time since explosion^ (yr) 


~ 20.1 


~ 16 


< 2 


~ 8.6 


(i/^SN 1993j)6 cm.peak 


~ 1.6 


~ 13 


~ 2 


1 


Radio brightness structure 


Shell (likely) 


Distorted shell 


? 


~Smooth shell 


Af/10- 5 M Q yr' 


~ (12- 16) 


~ 20 


? 


~ 5 


Deceleration parameter (m) 


~ 0.62 


-0.90 


? 


-0.82 


fb,eak (years) 


6 + 2 


Not yet 


~? 


-0.5 


Asymmetric expansion? 


No 


Yes 


? 


No(< 5%) 


Circumstellar medium 


? 


Clumpy 


? 


Approx. smooth 


M^ept/Mo 


~ 1.6 


~ 2.2 


~? 


-0.4 


M env /M 


~ 0.9 


~ 12 


~? 


-0.6 


Explosion scenario 


Binary 


Single 


Binary? 


Binary 


Magnetic field amplification 


Turbulent 


Turbulent 


Turbulent 


Turbulent 



1 From Marcaide et al. 1 2002 1. 2 From Perez-Torres et al. 1 2002 1. 3 From Perez-Torres et al. 1 2004 1. 4 From Marcaide et al. 119951 . Marcaide 
et al. {1995:1, Marcaide et al. f\WT\ . Perez-Torres et al. <2001l2002t . 5 Time since explosion at which we are discussing the property of the 
supernova. 



(optically thin shell), up to 14900+3000 km s -1 (optically thin 
sphere). For our second epoch (f = 583 d ), we obtain veloci- 
ties between 7200+700 km s _1 (shell) and 8600+ 1800 km s _1 
(sphere). For any given model, the expansion speeds seem to 
decrease, but this decrement is only marginal, so the expansion 
is compatible with a non-decelerated ejecta. Thus, future VLBI 
observations are necessary to confirm, or reject, the validity of 
the deceleration. 

We can estimate a minimum total energy and a minimum 
magnetic field for SN 200 lgd, as we have done previously for 
SN 1986J. From our VLA observations on 8 April 2003, we 
determined an spectral index of a = -1.05 + 0.08 (S v oc v") 
between 1 .4 and 43 GHz, and S 8 .4gh z = 1 .02 + 0.02. With those 
values, we obtain L v = (6.0 + 0.1) x 10 36 Dj 3 erg s~'between 
1.4 and 43 GHz. As R, we take the linear size corresponding 
to half the angular size of SN 200 lgd measured in the 8 April 
2003 epoch, 6 = (190 + 40)//as, which translates into to a linear 
size of R = D ■ G = (3.8 ± 0.8) x 10 16 cm. With these values, 
we get 

E mm » 1.6 X 10 46 (1 + kfl 1 $g 0% 7 O D\] n erg 

Bmin * 43(1 +k) 2 ' 1 <f> f6-^ 1 D^ n mG 

Since 1 < k < 2000, E min is in the range (2.4 X 10 46 - 
1.3 x 10 48 )erg, and the equipartition magnetic field in the 
range (52-376) mG. Thus, for a typical energy of 10 51 erg for 
the explosion of SN 200 lgd, the fraction of energy neces- 
sary to power the radio emission of SN 200 lgd is quite mod- 
est. The upper range of the equipartition magnetic field is in 
marginal agreement with that obtained for SN 1993 J at the 
same radius (B « 1.8G,[T33S| B * 0.5 G, USB. We note that 
our inferred equipartition magnetic field cannot be explained 
by compression of the pre-existing field in the circumstellar 
wind. Following the same reasoning made for SN 1986J, we 
obtain that B csm < (M v vl ,) 1/2 r~ l « 8 (M_ 4 vi ) 1/2 r^ 1 mG For 
r — 3.8 x 10 16 cm and a mass loss rate similar to that of 
SN 1993J (M = 5 x 10~ 5 I^yr 1 ), we obtain B csm < 1.5 mG, 
which is a factor about 35 to 250 times smaller than the equipar- 



tition field. Thus, it seems that to explain the radio emission 
from SN 200 lgd, we also need to invoke amplification mech- 
anisms other than compression of the circumstellar magnetic 
field, e.g., turbulent amplification. 

4. Summary 

We have presented results from VLBI observations of two Type 
II supernovae, SN 1986J and SN 200 lgd. Those results are 
summarized in Table 1, where we have also included for com- 
parison results obtained from radio observations on two other 
Type II supernovae, SN 1979C and SN 1993 J. 

VLBI observations of Type II supernovae suggest that if 
particles and fields are not far from equipartition, then it seems 
clear that amplification mechanisms other than compression 
from the existing circumstellar magnetic field must be acting, 
in order to explain the observed synchrotron radio emission. 
Turbulent amplification seems to be the most promising mech- 
anism. 

We would like to end this contribution by stressing that 
very-long-baseline interferometry (VLBI) observations of ra- 
dio supernovae are a powerful tool to probe the circumstellar 
interaction that takes place after a supernova explodes. We re- 
call the VLBI community that we have had to wait 1 1 years 
to have a really "SN1993J"-like event, and even so at a dis- 
tance four times greater. Since those supernovae belong to the 
rare group of radio supernovae that permit their follow-up with 
VLBI, we cannot miss the few opportunities we have of observ- 
ing and monitoring them. Unexpected results, like the discov- 
ery of a pulsar nebula or the accretion of material onto a black 
hole, might be waiting for us to unveil them, sometimes many 
years after the supernova explosion (e.g. Bietenholz, Bartel & 
Rupen 120041 . 
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